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ABSTRACT 


Halon 1301 is a halocarbon fire extinguishing agent (CBrF 3 ) used by the 
U.S, Army for vehicle fire suppression. Halon 1301 is discharged under nitrogen 
pressure, and the Halon-nitrogen mixture is a two-phase, two-component mixture 
that obeys compressible fluid laws and exhibits choking effects. A computer 
model was developed for analyzing the discharge of Halon and nitrogen from a 
storage bottle through pipes and nozzles. The model agrees well with data from 
Halon 1301 discharge tests. According to the model the discharge time depends 
mainly on nozzle area and pipe volume, for given bottle initial conditions. 
Graphs were developed for estimating discharge times. In tests of dispersion 
nozzles, conventional orifice-type nozzles gave only narrow-angle coverage, but 
a nozzle employing multiple concentric converging-diverging nozzles was 
developed which gave hemispherical coverage. 
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I. INTRODUCTION 


, ... - M s A rmy to fire extinguisher agents 

"Hal on” is the name given by 

, R f u "Halon 1301" is bromotr if luoromethane, CBrF 3 . 
containing halogens (Ref. 1). Halon 

arhnn three fluorines, no chlorine, and one 
The numerals "1301" stand for one carbon, 

bromine. A Halon 1301 concentration of only 7 percent by volume xn - 
extinguishes Tines by a cognation of cooling and intense 
chemical reaction chain of fuel and oxygen (Refs. 2-5). 


Halon 1301 is a clear liquid with a vapor pressure of 1.« MPa (200 psi ) 

room temperature. To speed expulsion, Halon is stored under nitrogen pressure 

„ c p MPa (750 psia) at 2U o 
The A r my has adopted a standard storage pressure of 5.2 

(Ref. 6). 


Figure 1 is a photograph of Halon 1301 discharging 

. wap. bottle valve restricting the flow as in 
The liquid discharge time, with only the bottle 

. . „ _ f nitroeen and Halon vapor takes 
this test, is 0.17 s. The subsequent venting 

nnff of white vapor. High speed 
another 0.2 s. To the eye the discharge is P 

movies show an opaque white plume that steadily diminishes 

change in density. In particular, there is no change in appearance, 

front lighting or back lighting, when the liquid flow stops and vapor venting 

begins, because the vapor condenses and continues to form a vapor cloud. 


..„,i -t of only a bottle and a nozzle, 
A Halon fire extinguisher system may consist 

or 3 complex arrangement of pipes and noxxles. m any oase. the design problem 
u to determine the bottle, noxxie. and pipe sixes that will deliver the desrred 
.mount of Halon at each discharge point in the desired time. A necessary design 
tool for this task is a flow model for predicting Halon flow rate versus t me as 









* K 



Figure 1. Halon 1301 discharge 


a function of bottle* pipe, and nozzle 


sizes. Such a flow model was the primary 
goal of this project. In addition, the flow distribution from dispersion 
nozzles was investigated. Tests were also conducted on valves, burst discs, 

flow deflectors, safety plugs and pressure gages; those te.ts were reported in 
Ref. 7. 


One benefit of a theoretical flow model is to provide a detailed 
understanding of the processes occurring during Halcn flow. Section II 
describes how nitrogen-pressurized Halon expands out of a bottle and through 
pipes and nozzles, based on the model. Section II also shows the verification 
of the model by comparison of theoretical and experimental pressure-time curves. 

Section III shows how to calculate Halon flow rate using the computer 

program and how to estimate discharge time using graphs calculated by the 
program. 


Section IV describes tests of Halon dispersion nozzles. Appendix A gives 
the derivation of the theoretical model. Appendix B presents equations for the 
properties of Halon and nitrogen mixtures. Appendix C describes work done 
toward modifying a steam blowdown program as an alternate approach to 

calculating Halon flow. Appendix D summarizes the bottle discharge tests made 
to provide a data base. 
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II. Nil ROGEN-PRESSURJZED HALON FLOW 


The purpose of this Section is to describe the behavior of Halon flow 
expanding fro. a bottle and flowing through pipes and orifices, based on the 
theoretical model. Comparisons between theoretical and measured pressure-time 
curves will be used to validate the model. 

A. Bottle Discharge Tests 

Figure 2 shows a typical pressure-time curve for a Halon discharge test 
(Test 116) and the apparatus used. The bottle had an electric pilot-operated 
valve. Nozzles or pipes could be connected to the valve exit. A pressure 
transducer and a thermistor temperature probe were mounted at the top of the 
bottle, and another temperature probe was mounted at the bottom. 


Before each run the bottle was loaded with the desired mass of liquid Halon 
(typically filling 60 percent of the bottu) a „ d pressuri2M wlth nltroge „_ ^ 

bottle was rocked and nitrogen added until the pressure stayed constant at the 
desired level, showing that the dissolved nitrogen concentration (about 2.3 
percent by mass) was in equilibrium with the nitrogen gas. 


The recorders were started, the valve was actuated 

temperature traces were recorded. The pressure dropped to 

ranging from a fraction of a second to several seconds, 
size. 


and the pressure and 
atmospheric in times 
depending on nozzle 


In most tests the 
key events in the run. 
The first slope change 


pressure-time curve had visible slope changes marking two 
The curve in Fig. 2 shows these slope changes clearly, 
is a pressure recovery from 3.9 MPa to 4.1 MPa at 0.2 s. 






The second change is an increase in pressure decay rate at 0.85 s. The pressure 
recovery at 0.2 s is due to a sudden release of dissolved nitrogen. The slope 
change at 0.85 s marks the point at which the bottle runs out of liquid and 
venting of nitrogen gas and Halon vapor begins. The identification of these 
slope changes with nitrogen release and liquid runout is based entirely on the 
theoretical model; there was no experimental observation of nitrogen release or 
liquid runout. 

The measured temperatures are shown in Fig. 3. The temperature measured by 
the top probe is considered to be the gas temperature, and the temperature 
measured by the bottom probe is considered to be the liquid temperature. The 
measured temperatures drop gradually until liquid runout at 0.85 s. Then the 
temperatures drop more rapidly until limited by heat input from the warm bottle 
walls . 

B. Bottle Discharge Predictions 

The theoretical flow model predicts the pressure decay curve shown in Fig. 
4 for Test 146. Numbers on the curve identify several events predicted by the 
model. During the first part of the run, from Point 1 to Point 2, nitrogen 
stays dissolved in the Halon as a non-equilibrium supersaturated solution, and 
the Halon remains a clear liquid. At Point 2 the pressure has dropped 
sufficiently to allow nitrogen bubbles to form. Once the bubbles form, all of 
the excess dissolved nitrogen quickly comes out of solution into the bubbles. 
The bubbles raise the Halon liquid level by about 25 percent and recompress the 
ullage gas to the pressure at Point 3. 


6 





PRESSURE. MPa 






The last liquid leaves the bottle at Point 4. The rate of pressure decay 
then increases because the gas is less dense than the liquid and allows a higher 
volume flow rate through the nozzle. 

"Outage fraction" is the mass fraction of the original Halon fill that has 
been discharged from the bottle at any given time. Figure 5 shows the 
theoretical increase of outage fraction with time. At liquid runout the outage 
fraction is 0.93; only 7 percent of the Halon remains in the bottle. 

C. Nitrogen Release Pressure 

The pressure at which the nitrogen starts coming out of solution, the 
"nitrogen release pressure" (identified by Point 2 in Fig. 4), depends on the 
Halon surface tension. As the pressure drops from Point 1 to Point 2, 
microscopic nitrogen bubbles are continually forming throughout the liquid due 
to random agglomeration of nitrogen molecules. The pressure of the nitrogen in 
these nucleation bubbles is equal to the initial nitrogen pressure, because the 
dissolved nitrogen concentration has not changed. However, the bubbles do not 
grow when the bottle pressure drops, because the surface tension of the Halon 
produces an additional inward pressure. Bubble growth only starts when the 
bottle pressure falls below the nucleation bubble pressure by an amount equal to 
the surface-tension pressure. Once this bottle pressure is reached the bubbles 
start to grow, and the surface-tension pressure on the bubbles decreases due to 
the increasing radius of curvature. The bubble growth accelerates, and the 
excess nitrogen in solution is quickly released. 

The difference between nucleation bubble pressure and bottle pressure at 
the start of nitrogen release (Point 2) is plotted as a function of Halon 
temperature for several runs in Fig. 6. As shown by the inset pressure-time 



TIME, s 


Figure 


Theoretical outage fraction 



RELEASE PRESSURE DIFFERENCE & 





curve in Fig. 6, the nucleation bubble 


pressure, which is the sum of the Halon 



vapor pressure and the initial nitrogen pressure, decreases slightly during 

expansion from Point 1 to Point 2 because the Halon vapor pressure decreases. 
The pressure difference 4p rel at the start of „ itrogen relea3e de< , rea3e3 ulth 

increased Halon temperature, because the surface tension decreases. At 20°C the 
release pressure difference is about 1.5 MPa. At the critical temperature of 
Halon 1301, 67 C, the surface tension is zero and the observed pressure 

difference Ap rel a t nitrogen release is also zero. There is considerable 
scatter in the data; the pressure difference required for nitrogen release is 
three times as great in some runs as in others at the same temperature. 


The pressure difference due to the surface tension of the liquid 
surrounding a bubble is 4<x/ D, where <x is the surface tension and D is the bubble 
diameter. As shown in Fig. 6, a nucleation bubble diameter of 15 nm gives a 
pressure difference that agrees with the lowest release pressure differences 
observed. That value of nucleation bubble diameter is used in the model, 
because it also gives the best agreement with the data between 20°C and 30°C 
where most of the tests were run. 


D. Theory and Data Comparison 

For comparing theoretical and experimental pressure-time curves it is 
necessary to match the starting times. The procedure adopted is to match the 
curves a short time after the valve-opening transients are over, about halfway 
between Point 1 and Point 2. 




Figure 7 compares the theoretical and experimental pressure-time curves for 
Test 1H 6. The curves agree well up to the point of liquid runout, Point 4. 
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During subsequent vapor venting the measured pressure falls more slowly than the 
theoretical pressure. The slower decay could be due to liquid draining from the 
bottle walls and raising the density of the nozzle flow. 

The theoretical and experimental temperatures are compared in Fig, 8. In 
the model the gas and liquid temperatures are equal and there is no heat 
transfer from the bottle. The theoretical temperature falls more rapidly than 
the measured temperatures especially auring vapor discharge. During vapor 
venting the theoretical temperature drops to -BC^C whereas the measured 
temperatures remain above -20°C. The warmer measured temperatures are 
attributed to heat transfer into the gas and into the temperature probes from 
the bottle walls. 

Figure 9 compares the theoretical and experimental pressure traces for Test 
102, a test with no nozzle, using only the valve to restrict the flow. The 
valve flow area cannot be measured directly because of the complicated shape of 
the valve flow passage. Instead, the flow area used in Fig. 9 is the one that 
gives the best agreement between the theoretical and experimental pressure 
curves, 500 mm 2 . The flow area measured in water flow tests was less, 400 mm 2 . 
Test 102 can be considered a calibration of the effective flow area of the valve 
with Halon. The question of determining flow areas for use in the theoretical 
model will be discussed further in Section III. 

Because of the substantial agreement between the theoretical and measured 
pressure-time curves for Test 146, and other tests to be discussed (especially 
up to the time of liquid runout), it is concluded that all other quantities 
calculated by the model are also correct. It is unlikely that quantities such 


TIME, s 


Figure 8 


Comparison between theoretical and experimental temperature-time 
curves 






PRESSURE. MPa 


Figure 9 



Comparison, with valve only, 
pressures 


between theoretical and experimental 




ure, and nitrogen concentration could be 


as void fraction, quality, temperat 
significantly in error and the pressure-time curve still be correct. After 
liquid runout, the model is less accurate, as evidenced by the slower pressure 
decay and lower temperatures than measured. 


E. Bottle Expansion Process 

The theoretical model calculates many details of the Halon discharge 
process in addition to the pressure, temperature, and outage values that have 
been presented. A study of these other quantities provides useful insight into 
the nature of the flow. The processes occurring inside the bottle, as 
calculated by the theoretical model, will be considered first. 


Figure 10 presents sketches of the interior of the Halon bottle at four key 
times during discharge, with corresponding points marked on pressure and outage 
curves. Initially, at Point 1, the bottle is about 60 percent full of liquid. 
The ullage volume above the liquid contains nitrogen gas and Halon vapor. The 
region below the liquid surface, which will be re: -rred to as the "liquid 
layer," consists of liquid Halon and dissolved nitrogen. The dissolved nitrogen 
initially has a concentration of about 2.3 percent by mass. When the valve is 
opened, fluid flows from the liquid layer through the valve and nozzle to 
atmosphere. Nitrogen bubbles form somewhere between the valve entrance and the 

nozzle exit as the pressure in the flow stream falls below the nitrogen release 
pressure. 


As the bottle empties, Halon evaporates and supplies vapor to the growing 
ullage volume. The nitrogen dissolved in the evaporated Halon also contributes 
a small amount of additional nitrogen to the ullage. In addition, the dissolved 
nitrogen that is very close to the surface (within a few mm) probably comes out 
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Figure 10. Bottle expansion steps 


of solution and enters the ullage, but this effect is ignored in the model. The 
bulk of the dissolved nitrogen is too far below the surface to diffuse to the 
surface and is only released when bubbles form throughout the liquid at the 
nitrogen release pressure. Sketch 2 in Fig. 10 shows the conditions just before 
nitrogen release and Sketch 3 shows the conditions just after equilibrium 
concentration is reached. 

During bubble formation the liquid layer expands by about 25 percent, 
compressing the ullage gas and raising the pressure to Point 3. The liquid 
layer is now a two-phase mixture of (1) a liquid phase consisting of Halon and 
dissolved nitrogen and (2) a gas phase (distributed as bubbles) consisting of 
nitrogen gas and Halon vapor. 

The liquid layer continues to discharge from the bottle until all of the 
liquid has left the bottle at Point 4. At that time the ullage gas starts to 
flow out of the bottle. As the gas in the ullage cools, the Halon vapor starts 
to condense, and the flow leaving the bottle becomes an increasingly wet mixture 
of nitrogen gas, Halon vapor, and Halon liquid. 

The conditions inside the bottle depend only on the amount of Halon that 
has been discharged, and not on time, under the assumption of adiabatic flow. 
Figure 11 shows how the pressure and temperature decrease as the outage fraction 
increases, for typical initial conditions. ~he pressure drops to 3.8 MPa at an 
outage of 0.30 before the first nitrogen bubbles form (based on a 15-nm 
nucleation bubble diameter). The temperature drops to 15°C at the same time. 
As shown in Fig. 12, the mass of liquid in the bottle (Halon plus dissolved 
nitrogen) has dropped to 2.1 kg at this point. The mass of gas in the ullage 
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has risen from an initial 0.22 kg to 0.29 kg because of Halon evaporation. 

When the nitrogen bubbles form, the pressure increases by 0.4 MPa, the 
temperature increases by 1°C, the liquid mass decreases by 0.1 kg, and the 
ullage gas mass decreases by 0.03 kg (due to Halon vapor condensation). The 
mass of bubbles produced is 0.05 kg. As shown in Fig. 1 3» the exit flow becomes 
a two-phase mixture having 0.02 quality and 0.2 void fraction. 

The mass of gas in the liquid layer reaches a maximum of 0.06 kg at 0.49 
outage. The mass of gas in the ullage reaches a maximum of 0.30 kg at 0.86 
outage. The last liquid ieaves the bottle at an outage of 0.90; 10 percent of 
the Halon remains as vapor. The quality of the last of the liquid layer leaving 
is 0.07 and the void fraction is 0.60. 

After the last liquid leaves the bottle the ullage gas follows a pressure 
and temperature path that causes condensation of Halon vapor. The quality of 
the exit flow starts at 1.0 at liquid runout and drops to 0.75 as the pressure 
drops to atmospheric. The void fraction only drops to 0.998. The final outage 
fraction is 0.995; 0.5 percent of the Halon remains in the bottle at atmospheric 
pressure. 

F. Nozzle Flow 

The flow rate out of the bottle depends on the response of the external 
flow components to the bottle conditions. An important case is a nozzle 
connected directly to the valve. The nozzle may have multiple passages to 
disperse the Halon, but the net effect is an area restriction through which the 
Halon must pass. As the area decreases in the flow direction, the stream 




pressure drops, but the pressure oannot reach atmospheric in a converging 
passage because the flow is compressible and sonic conditions are reached at 
about half the nozzle inlet pressure. Sonic flow is reached at the minimum-area 
point, or "throat" in the nozzle. The nozzle can either be cut off at that 
point for unguided fluid expansion to atmospheric pressure (as in Fig. 1 . where 

the valve acts as a converging nozzle) or the nozzle can have a diverging 
section for guided expansion. 


Fisure Id shows the calculated flow conditions in a typical converging- 
diverging nozzle (or in a converging nozzle cut off at the throat). The inlet 
conditions are the bottle conditions of Fig. „ at 0.5 outage fraction: 3.6 MPa. 
12°C. 0.04 quality, and 0.34 void fraction. The pressure decreases to 2.4 MPa 
at the throat. The quality increases to 0.09 and the void fraction to 0.06. 
The velocity at the throat is 56 m/s. The throat temperature is 7«c. 

As the flow proceeds through the diverging section (or to atmosphere in an 
unguided expansion) the quality, void fraction and velocity continue to 
increase. When atmospheric pressure is reached (at a flow area 11 times the 
throat area) the quality is 0.47. the void fraction is 0.996. and the velocity 
is 200 m/s. If the nozzle is cut off at the throat, the same conditions would 
be reached except that the gas phase may expand laterally away from the liquid, 
leaving a slow liquid core surrounded by a higher-velocity gas sheath. 

The flow rate set by the 314 mm2 throat area at 3.6 MPa nozzle inlet 
pressure is 10.4 kg/s. The flow rates at other times during the bottle blowdown 
are plotted as a function of time in Fig. 15 . The initial flow rate is 15 kg/s. 
me flow rate increases when the pressure rises at the time of nitrogen release. 
The flow rate drops abruptly when the liquid runs out and gas venting begins. 














G* Pipe Flow 


Another Important type „ f oomponeot ^ p plpe ^ ^ ^ 

to atmosphere or to nozzles at a distance. The leading edge of the Halon 
tracing down the pipe forms a "liquid front" oonslsting of a two-phase Halon- 
nitrogen mixture traveling at sorie velocity. If the pipe leads to atmosphere 


Figure ,6 shows the theoretical variation of pressure along a 314 mm? pipe 

at the same inlet conditions as in Fig 14 T h« „ al , * . 

&• . The calculations apply to the 

instant the liquid front has traveled 2 4m d nun a i 

d.n m down a long pipe or when steady flow 

has been established at the end of a 2.4-m pipe. 


The pressure falls to 3.3 MPa at the 
liquid front or pipe exit. The velocity is 
at the exit. The main acceleration takes 
liquid front or pipe exit. 


pipe inlet and to 1.7 MPa at the 
25 m/s at the pipe inlet and 59 m/s 
place as the flow approaches the 


The flow rate set 
reduction from the 10.4 
mm 2 , 2.4-m pipe fills, 
kg/s. 


by choked flow In the pipe is 6.9 kg/s. a 34 percent 
Ag/s with a nozzle of the same area. Thus, as a 31. 
the flow rate starts at 10.4 kg/s and decreases to 6.S 


H. Pipe Pressurization 

If there 13 3 002216 ° r 0th - -■‘'‘'“on at the end of a pipe, the pipe 
exit pressure starts to rise when the liquid front reaches the end. The 

pressure rises until the flow rates into and out of the pipe are balanced at a 
lower flow rate. The theoretical pressures along the pipe at various stages 
during the filling and pressurization process are illustrated in Fig. 17, for 

2 r 



VELOCITY V, m/s 
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Halon discharging from a 3750 cm3 bottle into a 2550 cm3 pipe with a 360 mm^ 
nozzle at the end. The initial bottle pressure is 5.2 MPa. When the valve is 
opened the Halon quickly reaches the pipe inlet and establishes sonic flow at 
that point; the pressure profile is shown by Curve 1. The pressure at the valve 
exit at that time is 4.1 MPa and the pressure at the pipe inlet is 2.8 MPa. The 
flow rate is 24 kg/s. 


The liquid front proceeds down the pipe at 60 m/s and reaches the end in 62 
ms. As shown by Curve 2, the bottle pressure has dropped to 3.6 MPa, but the 
pipe inlet pressure has increased to 3.3 MPa. The pressure at the liquid front 
(on the liquid side) is 1.6 MPa. The flow rate is 14 kg/s. The mass of fluid 
in the pipe is 1.5 kg, 47 percent of the Halon fill. 

The Halon then encounters a nozzle with a throat area about half the pipe 
area. The flow decreases and the pipe exit pressure increases. The exit 
pressure reaches a peak 31 ms after the arrival of the liquid front (Curve 3). 
The bottle pressure at that time is 3.2 MPa, the pipe exit pressure is 2.7 MPa, 
and the flow rate is 9 kg/s. The amount of Halon stored in the pipe is 1.9 kg, 
60 percent of the Halon bottle fill. About 6 percent of the Halon has been 
discharged from the nozzle by this time. 

I. Fipe and Nozzle Discharge Tests 

To verify the theory, flow tests were made with a 670-mm I. 2 , 3.8-m long pipe 
and various nozzles. Figure 18 is a sketch of the apparatus. The theoretical 
and experimental pressure curves, for both the bottle pressure and the pipe exit 
pressure, are compared in Fig. 19 for Test 175, in which a 360 mm 2 nozzle was 
used. The agreement is good. Several events predicted by the model are marked 
along the curves. At Point 1 the liquid front reaches the end of the pipe and 
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the pipe exit pressure starts to rise. In the test the exit pressure rises 
earlier due to air pushed ahead of the Halon. 

Point 2 is the instant of peak pipe exit pressure. At Point 3 the last 
liquid has left the bottle, and vapor is following the liquid down the pipe. At 
Point 4 the last liquid has left the pipe. The abrupt change in slope at liquid 
runout predicted by the theory is not observed in the test. 

The theoretical and experimental temperatures at the end of the pipe are 
compared in Fig. 20. The model predicts an initial temperature reading of 0 °C 
when the liquid front reaches the end of the pipe, followed by an increase to 
12°C as the pipe pressurizes. The measurements, however, show a 50°C spike due 
to the pulse of air preceding the Halon. The measured temperatures during gas 
venting (after Point 4) do not drop as low as predicted by the model, presumably 
because of heat transfer from the bottle and pipe. 

Theoretical and experimental pressures for Test 174, with a smaller nozzle, 
are compared in Fig. 21(a; for the bottle pressure and Fig. 21(b) for the pipe 
exit pressure. The characteristic differences are again evident: no sharp break 
in slope at liquid runout in the test and a longer venting time. Also, the 
measurements show oscillations in the pressure after the liquid reaches the end 
of the pipe; in fact, the pipe exit pressure is alternately higher and lower 
than the bottle pressure. 

Fig. 22 compares the theoretical and experimental pressures for a test with 
no exit nozzle. The flow expands directly from the pipe exit. Events 1 and 2 
(the liquid front reaching the exit and the pipe reaching peak pressure) are 













simultaneous. The agreement between theory and experiment is still good in this 
e.-treme case of a large pressure drop in the pipe. 

J. Multibranch Flow Systems 

A Halon distribution system may consist of several pipes and nozzles 
delivering Halon to multiple locations. The computer model was developed to the 
point of handling a single pipe leading to a manifold, with any number of pipes 
leading from the manifold to discharge nozzles. 

Two double-branch distribution systems were assembled to test the model. 
Figure 23 shows System 1 which had a 3.73-m pipe leading from the bottle to a 
tee and two pipes of different lengths leading to nozzles of different 
diameters. The inlets to the pipes, tee, and nozzles had rounded entrances to 
avoid any uncertainty about flow area. 

Figure 24 presents the theoretical pressure-time curves for the bottle, 
tee, and nozzles, and Fig. 25 shows theoretical Halon mass discharged versus 
time from each nozzle. As shown in Fig. 24, the liquid front reaches the 
nozzles very quickly. Traveling at about 60 m/s, the Halon reaches the most 
distant nozzle, 8 m from the bottle, in only 0.13 s. The piping system reaches 
peak pressure only 0.2 s after valve opening. Once nozzle flow begins, there 
are only small pressure drops in the pipes: 3 percent from bottle to tee, 7 

percent from the tee to the small nozzle, and 13 percent from the tee to the 
large nozzle. 

The bottle runs out of liquid at 1.0 s and the main pipe at 2.1 s. The 
short branch runs out of liquid at 2.3 s, clearing a vapor path from bottle to 
atmosphere and starting a more rapid pressure decay. Because of the diminishing 
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pressure, the longer branch does not run out of liquid until 2.8 s. 

The amount of Halon theoretically discharged from each nozzle (Fig. 25) is 
proportional to the nozzle area. At the time of liquid runout in the short 
branch, 89 percent of the Halon has been discharged. At the time of liquid 
runout in the long branch, 96 percent of the Halon has been discharged. 

Figures 26(a)-(d) compare the theoretical pressure-time curves with the 
measured pressure-time curves. The general agreement is good. Specific 
differences are that the theoretical bottle pressure does not show the observed 
undershoot when the system fills, and the theoretical pressures fall slightly 
below the data for most of the run, especially during vapor venting. 

The second multibranch system tested is sketched in Fig. 27. This system 
was designed to be an extreme example of a large length and high volume. The 
distance from bottle to nozzles is 9.8 m, and the piping volume is 6 percent 
greater than the bottle volume. The multibranch option of the computer program 
would not run with such a large volume without further programming effort. 
Instead, the system was analyzed as an equivalent single-branch system using the 
rules for equivalent pipe and nozzle sizes presented in the next Section. 

Figure 28 compares the theoretical and experimental pressure-time curves 
for this multibranch system 2. Because of the large pipe volume the bottle runs 
out of liquid before the pipes are fully pressurized. Nozzle flow does not 
start until all of the Halon has been transferred to the pipes and the bottle 
pressure has dropped to 2.1 MPa. The agreement between the theoretical and 
experimental pressures is reasonably good, with the theoretical pressures again 
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Figure 27. Two-branch distribution system 2 
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being somewhat lower than experimental. 



K. Summary 

The results presented in this Section have shown that the theoretical model 
predicts the pressure-time behavior of Halon discharge reasonably well. Other 
behavior predicted by the model, such as mass discharged versus time, is, 
therefore, probably also correct. Based on the model, a clearer picture of 
Halon discharge behavior emerges. After valve opening, the Halon travels 
rapidly to the discharge nozzles, and rapidly pressurizes the piping system. 
These steps are completed before any significant amount of Halon has been 
discharged and in a time that is short compared with the total discharge time. 
During discharge, the flow rates through the individual nozzles are divided in 
proportion to nozzle throat areas. The discharge time with a pipe between the 
bottle and nozzles is increased, compared with the discharge time for u nozzle 
only, by both pipe friction and by the loss of pressure due to filling the 

pipes . 



III. FLOW CALCULATION METHODS 


Two methods of calculating Halon and nitrogen flow were developed during 
this project. One method was a modification of a steam-and-water blowdown 
program developed at the Los Alamos National Laboratory (LAND, which was run on 
the LANL computer. Results of that approach are discussed in Appendix C. The 
other method was a program called HFLOW, written specifically for Halon and 
nitrogen discharge, and run on the JPL computer. The two methods agreed, but 
the HFLOW program was more convenient to use and was the model used for the 
theoretical calculations that have been presented in this report. 

.he multibranch version of HFLOW is limited in the piping geometries it 
will handle. However, multibranch systems can be reduced to single-branch 
systems by rules that will be presented, and the program runs well for single 
branches. For estimation of discharge time it is possible to use the computer 
program to prepare graphs of discharge time as a function of pipe length for 
various pipe and nozzle diameters. Such a graph can be prepared for each bottle 
condition of interest, and an example will be presented. Further simplification 
is possible if the flow system uses optimum pipe sizes. A time-estimation graph 
and formula will be presented for that case. 

A. Program HFLOW 

To run HFLOW for a single-branch system the following inputs are required: 
bottle volume, Halon mass, initial pressure, initial temperature, valve flow 
area, pipe length, pipe area, and nozzle area. The program prints various 
quantities at each of 100 pressure steps, including time, flow rate, mass 
discharged, mass stored in the pipe, and temperatures and pressures in the 
bottle and at the pipe exit. 







Determining the flow areas to be used in the program may be difficult for 
some types of nozzles and certainly for the bottle valve. The flow area can be 
known accurately only for orifices having a well-rounded entrance, with a radius 
equal to at least half the orifice diameter, in which case the flow area can be 
assumed to be equal to the physical area. 

If the orifice has a sharp entrance or only a small radius, the effective 
flow area can be as little as 60 percent of the physical area. In such cases 
the flow area (or the discharge coefficient, which is the ratio between the 
effective flow area and physical area) can be determined from flow tests with 
the type of fluid to be used. Discharge coefficients for some orifice shapes 
can be found in hydraulics handbooks. 

Flow areas determined from water or gas flow tests or from handbook 
coefficients can be in error by 10 to 20 percent for two-phase flow, and the 
calculations using HFLOW will reflect that error. However, for flow components 
upstream of the final discharge nozzle, such as the valve, the effect of flow- 
area errors is usually small. If the discharge nozzle has a flow area less than 
70 percent of the upstream restrictions, 90 percent or more of the pressure drop 
from the bottle to atmosphere will take place across tne discharge nozzle, 
making it only necessary to provide the program with an accurate flow area for 
the discharge nozzle. 

If the final discharge is through a complex restriction, such as the valve 



itself, then either a water calibration can be used, with uncertain accuracy for 
two-pha 3 e flow, or a Halon discharge test can be used a3 a calibration. As 



shown earlier (Fig* 9), the best-fitting valve flow area for a Halon discharge 
test was 500 mm2 whereas the valve flow area from a water calibration was 
400 mm^. 

The equations used in program HFLOW are derived in Appendix A. The overall 
organization of the program is summarized in Fig. 29. The program has two 
sections. The first section calculates bottle conditions as a function of 
bottle pressure. The second section calculates flow rates and times. 

The bottle section has a main program that calls subroutines for fluid 
properties, nitrogen release, and equation solving. The program proceeds 
stepwise through 100 pressure increments from the initial bottle pressure to 
atmospheric pressure. At each pressure step the temperature of the fluid in the 
bottle is guessed. The subroutines calculate the properties of the Halon and 
nitrogen mixture at that temperature and pressure. The main program calculates 
the amount of fluid that must be discharged in each step to reach the new 
pressure. The energy equation is then checked: change of internal energy of the 
bottle fluid equals enthalpy of the flow discharged. If the two quantities are 
not equal the subroutine SFABZ (subroutine to find a bounded zero) is called to 
furnish a new temperature guess. Iterations continue until the correct 
temperature is found. 

When the pressure reaches the nitrogen release pressure, based on 15 nm 
nucleation bubble diameter, the nitrogen release subroutine calculates the 
conditions after transition to equilibrium nitrogen concentration at constant 
internal energy. The pressure steps then continue. 

The result of the bottle section of the program is a table of the fluid 















mass increments discharged in each pressure step, and values of the pressure, 
temperature and quality of each mass increment discharged. These results are 
used by the flow section of the program. 



The flow section of the program has a main program th calls subroutines 
for fluid properties, component pressure drops, and equation solving. At each 
pressure step the main program guesses a flow rate. A downward pressure scan is 
made through each flow component (valve, pipe, and orifice) that is carrying 
flow at that time. At each pressure step the temperature of the fluid at that 
point is guessed, the property routines are called, the energy equation is 
checked, and SFABZ is called to furnish a new temperature guess to satisfy the 
energy equation. When the point in the flow that is at sonic conditions is 
reached the required flow area is calculated and compared with the actual flow 
area of the flow component that is limiting the flow at that time (the pipe 
area, if the liquid front has not reached the end of the pipe, and the nozzle 
area after that). If the required and actual areas are different, subroutine 
SFABZ supplies a new flow rate guess. 

Once the correct flow rate is found, and added to the rate of fluid mass 
change in the pipe, the time required to discharge the required flow increment 
from the bottle is calculated and added to the elapsed time. If the liquid 
front has not reached the end of t» e pipe, the front is then advanced by its 
travel distance in that time increment. 

The result of the flow section of the program is a printout of elapsed time 
at each step of bottle pressure, together with values of temperatures and 
pressures, fluid front positions and other quantities in the bottle and pipes. 



B. Reducing Multibranch to Single Branch Systems 

Most distribution systems have relatively small pressure drops from bottle 
to nozzles once nozzle flow begins. Therefore, the pressure in the bottle and 
pipes at the start of nozzle flow depends mainly on the total pipe volume and 
not on how the volume is distributed throughout the system. In addition, the 
flow rate out of the bottle is the sum of the nozzle flow rates, plus the rate 
of storage change in the pipe volume, and the bottle discharge flow rate does 
not depend greatly on nozzle location. 

Since the locations of the pipe volumes and nozzle areas are of secondary 
importance, it should be possible to calculate the bottle and nozzle discharge 
flow rates for a multibranch system by treating the multibranch system as a 
single branch system with (1) the single branch having the same volume as the 
multiple pipes and (2) the single nozzle having the same total area as the 
multiple nozzles. 

The nozzle areas to be used are effective areas reduced from the actual 
nozzle areas by factors that account for pressure drop in the pipes. Reduction 
factors for Halon flow at typical mean bottle conditions during discharge (3.0 
MPa, 10°C and 4 percent quality) are plotted in Fig. 30 as a function of 
length/ d ia m eter ratio for pipes having areas of 1.5, 2, 3» and 4 times the 
nozzle area. 

An example of reducing a multibranch system to a single branch system will 
be presented for the system of Fig. 23. For Branch 1 the pipe/nozzle area ratio 
is 1.8 and the pipe L/D is 150. From Fig. 30 the area reduction factor is 0.90, 
giving an effective nozzle area of 47 mm 2 . For Branch 2 the area ratio is 3.7, 
the L/D is 400, the reduction factor is 0,96, and the effective area is 24 nm^. 
The total effective area of the branches is thus 71 mm 2 . 




For the main branch the ratio of pipe area to effective nozzle area is 
412/71 = 5.8 and the L/D is 160. From Fig. 30 the reduction factor is close to 



1.0. The total volume of the pipes is 2100 cm3, 
length and 5.1 m length has the same volume. 


A single pipe of 412 mm 2 



' Thus, a single branch system equivalent to the two-branch system of Fig. 23 

< consists of a 412-mm 2 , 5.1-m long pipe and a 71 mm 2 nozzle. When these values 

| are used as inputs to program HFLOW the result is a single curve of mass 

; discharged versus time. The mass discharged from each nozzle in the multibranch 

system is obtained by multiplying the total mass discharged by the fraction of 

total nozzle area for each nozzle, namely 47/71 for Branch 1 and 24/71 for 

. Branch 2. The mass discharge curves obtained in this way are compared in Fig. 

31 with the curves from the multibranch run of the computer program. The masses 

[ discharged are about the same for the two program runs, and the discharge time 

1 for the single-branch approximation is about the same as for the earliest- 

» 

emptying branch of the two-branch system. Thus, the single-branch approximation 
( appears to be a valid approach to simplifying the computations. 




C. Discharge Time Graphs 

To provide a convenient way of estimating discharge times, the computer 
program can be run for a range of pipe sizes and orifice sizes and the results 
plotted in a graph for each bottle size and set of initial conditions. Figure 
32 presents such a graph for a 3750 cm3 bottle with 3.2 kg of Halon at 5.2 MPa 
initial pressure and 20°C initial temperature. Discharge time is plotted as a 
function of pipe length for four nozzle areas and several values of pipe/nozzle 
area ratio. As would be expected, the discharge times decrease as nozzle area 
increases, and the discharge times increase with pipe length. 
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The discharge time for the two-branch system of Fig. 23 is estimated by 
using the effective nozzle area of 71 mm 2 , the pipe area of 412 mm 2 , and the 
effective pipe length of 5.1 m. Fig. 3 2 shows that for a 75 mm 2 nozzle area and 
a pipe/nozzle area ratio of 5*8 the discharge time at 5.1 m length is about 
2.0 s. For small changes in nozzle area the discharge time is proportional to 
area; thus, the discharge time for a nozzle area of 71 mm 2 is 75/71 x 2.0 = 2.1, 
which is close to the value given by the computer program in Fig. 31. 

For the flow system of Fig. 27. the effective nozzle area for each branch, 
from Fig. 30, is 105 mm 2 . The total pipe volume is 3990 cm3, the pipe/nozzle 
area ratio is 2, and the effective single-pipe length is 9*7 m. From Fig. 32, 
interpolating for 210 mm 2 area and extrapolating to 9.7 m length, the discharge 
time is about 1.2 s, in agreement with Fig. 28. 

Thus, a graph of the type shown in Fig. 32 provides useful estimates of 
discharge time. The computer program could be used to generate a set of more 
detailed graphs, with additional values of nozzle area, for each initial bottle 
condition of interest. 

D, Optimum Pipe Diameters 

It can be seen from Fig. 32 that there is an optimum pipe area, giving 
minimum discharge time, for each nozzle area. At larger pipe areas than optimum 
the discharge time is increased because of the pressure loss in filling the 
large pipe volume. At smaller pipe areas than optimum the discharge time is 
increased because of the friction pressure drop in the pipe. 

Figure 33 illustrates the effect of varying pipe/nozzle area ratio. 
Pressure-time curves are plotted for five different pipe areas at a fixed nozzle 








area of 77 mm2. In curve 1 the pipe has the same area as the nozzle and the 
discharge time is 2.4 s. Increasing the pipe area to 1.3 times the nozzle area 
gives a higher flow rate and a lower discharge time of 1.9 s. With a pipe area 
2.6 times the nozzle area, the discharge time is further reduced to 1.6 s. 
However, at a pipe area of 5.2 times the nozzle area, the loss in pressure due 
to filling the pipe begins to dominate and the discharge time increases to 
1.8 s. At a pipe area of 8.7 times the nozzle area, the discharge time is 
further increased to 2.1 s. Thus, the optimum ratio of pipe area to nozzle area 
is about 2.6 at this pipe length and nozzle area. 


Figure 32 shows that the optimum pipe/nozzle area ratio decreases with 
increasing nozzle area. The optimum area ratio ranges from 3.0 at 75 mm2 nozzle 
area to 1.5 at 300 mm2 nozzle area. However, an area ratio of 2.0 gives nearly 
minimum discharge time at all nozzle areas. 


If a Halon distribution system is built with optimum pipe areas (twice the 
nozzle area is close enough), then the estimation of discharge time can be 
further simplified to a single curve for each nozzle area. Furthermore, the 
discharge times can be normalized by nozzle area to factor out the first-order 
effect of nozzle area. Figure 34 shows such a graph. The product of nozzle 
area and discharge time is plotted as a function of pipe length for four 
different nozzle areas. For the Fig. 33 conditions of 77 mm 2 nozzle area and 
3.8 m pipe length, for example, Fig. 34 shows that the product of nozzle area 
and discharge time is 128 mm 2 -s. Therefore, the discharge time is 128/77 = 
1.7 s, in agreement with the minimum discharge time from Fig. 33 . 
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E. Discharge Time Formula 

The curves in Fig. 3^ can be approximated by straight lines fitted by the 
.ollowing formula: 

Discharge Time = 95 ♦ (*<.6 + 0.047 A)L seconds 

A 

where A is the nozzle area in mm^ and L is the pipe length in m. Fe>~ a 
multibranch system, area A is the sum of the effective areas of the branches 
(calculated from Fig. 30) and length L is the length of pipe having the total 
volume of the system. For the system of Fig. 27, for example, which has close 
to optimum pipe sizes, A = 210 mm2 anc | l = 9.7 m. The formula gives a liquid 
discharge time of 1.1 s, 15 percent less than the value given by the computer 
program in Fig. 28. 

F. Comparison With NFPA Method 

The National Fire Protection Asoociation has presented a method of sizing 
pipes and nozzles for Halon 1301 fire extinguishing systems (Ref. 8), and 
manufacturers have published handbooks for applying the NFPA method at 
particular initial pressures (Ref. 9). 

The procedures of the NFPA method are, in effect, the same as the ones used 
here in converting multibranch systems to single-branch systems: the nozzle 
areas (corrected for pipe pressure drop) are added to give a single nozzle area, 
and the pipe volumes are added to give a single pipe volume. 

Reference 9 presents a sample calculation for a Halon system with two 
bottles feeding six nozzles. The system is sketched in Fig. 35. The aim of the 
calculation in Ref. 9 is to size the nozzles and pipes for a "10 second 


discharge time". 
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From Fig. 30 the branches have area reduction factors from 0.82 to 0.92. 



The total effective area of the six nozzles is 2340 mm 2 . The choice of 
equivalent pipe length and diameter is arbitrary so long as the volume equals 
the total pipe volume; 8 m length and 7330 mm 2 area are chosen. Figure 35 shows 
the pressure-time curves predicted by program HFL0W, The time required for 
discharge of the entire amount of Halon, including vapor, is 7.2 s, which is 
consistent with the 10 s design value. However, the liquid discharge is 
complete at 4.3 s, and 90 percent of the Halon has been discharged by that time. 
Thus, it may be possible to use smaller pipes and nozzles than calculated by the 
NFPA method for a given discharge time. 

G. Summary 

The computer program HFLOW provides accurate calculations of Halon 
discharge flow for single-branch systems. Multibranch systems can be handled as 
equivalent single-branch systems. The results of computer runs can be plotted 
for convenient estimation of discharge times. 

The optimum pipe/nozzle area ratios for minimum discharge times are 1.5 to 
3.0, depending on nozzle area, and an area ratio of 2.0 will give near-minimum 
discharge times in all cases. For such ratios, the time, in seconds, required 
for discharge of the liquid (about 90 percent of the total Halon) from a system 
with a 3750 cm3 bottle containing 3.2 kg of Halon at 5.2 MPa and 20°C is given 
by L95 + (4.6 ♦ 0.047 A)LJ/A where A is the nozzle area in mm 2 and L is the pipe 
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length in m. 


IV. HALON DISPERSION NOZZLES 


A. Liquid Atomization and Acceleration 


About half of the mass of Halon leaving a discharge nozzle is liquid. The 


atomization and acceleration of this liquid affects the way the Halon travels 


after leaving the nozzle. Fortunately, two-phase nozzle flow is well understood 


(Ref. 10). The liquid acceleration mechanism is slip between the gas and liquid 


which, at the same time, atomizes the liquid to small drops. 


Applying the theory of Ref. 10 to the Halon nozzle discussed previously 


(Fig. 14), it is possible to replace the isentropic velocity plotted in Fig. 14 


with the actual gas and liquid velocities and to calculate the drop size. 


Figure 36 shows the results. The liquid atomizes to droplets of 15 pm diameter 


by the time the flow reaches the throat. With such small drops the liquid is 


accelerated efficiently, and the liquid velocity stays at about 90 percent of 


the gas velocity. 


B. Orifice Nozzles 


Halon nozzles used in the past have had converging or straight passages 


only, with no diverging section. In such "orifice nozzles" the exit velocity is 


only about 30 percent of the velocity attainable at atmospheric pressure, as can 


be seen from Fig. 36. The Halon continues to accelerate outside the nozzle, but 


with uncontrolled spreading and with less effective liquid acceleration than 


would be provided by a diverging section. 


To disperse the Halon into the compartment to be flooded, orifice nozzles 
use slots or holes aimed in the desired direction. The nozzles may also have 
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deflectors downstream to re-direct the flow leaving the nozzle. Three such 
nozzles were tested, and they are shown in Fig. 37. The "beehive" nozzle had 
slots cut in one side of a cylinder and a hole in the end. The "splitter" 
nozzle had two bars in front of a straight orifice to divide and spread the 
flow. The "showerhead" nozzle had 31 holes drilled at angles of 15, 45, and 60 
degrees from center in an effort to provide hemispherical coverage. The total 
orifice area was 1370 mm^ for the beehive nozzle, 790 mm^ for the splitter 
nozzle and 157 mm^ for the showerhead nozzle. 

The three nozzles were tested with water and nitrogen mixtures for visual 
and photographic observation and with Halon using high-speed movies. Figure 
38(a) shows the beehive nozzle operating with water and nitrogen. The flow 
leaves as a flat sheet from the slots on the side and as a jet from the end. 

Figure 38(b) is a movie frame of the beehive nozzle in a Halon test, at a time 
8 ms after the start of flow. The view is toward the hole in the end of the 
nozzle, showing the main Halon sheet flat on, traveling to the right. The 
beehive nozzle is effective only in producing a planar flow directed to one side 
of the nozzle. 

Figure 39(a) shows the splitter nozzle operating with water and nitrogen. 

The deflector bars produce four jet lobes spreading at about 20 degrees from 
center. Figure 39(b) shows the splitter nozzle with Halon at 2 ms from the 
start of flow. Four diverging Halon lobes can be seen. The splitter nozzle is 
effective in producing a slightly divergent discharge flow. 

The showerhead nozzle shown operating with water and nitrogen flow in 
Fig. 40(a). The flow initially 3pread3 at 60 degrees but pulls inward to a 
cylindrical spray downstream. With Halon flow, Fig. 40(b), the showerhead 
nozzle has an even more pronounced pulling-in effect, yielding a jet that has no 
more spreading than the flow from the bottle valve alone, as previously shown in 
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Figure 38. Flow test of the beehive nozzle: (a) With water and nitrogen 



Figure 38. Flow test of the beehive nozzle: (b) With ilalon 
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39 . Flow test of the splitter nozzle: (a) With water and nitrogen 



Figure 39. Flow test of the splitter nozzle: (b) With Kalon 
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Flow test of the showerhead nozzle: 


(a) With water and nitrogen 


Figure 
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Fig. 1. 

The reason that the showerhead nozzle fails to produce a spreading flow is 
simply that continuity prevails. At atmospheric pressure the Halon has a 
certain velocity and density that defines a fixed flow area. The Halon cannot 
expand to a larger area without going to lower than atmospheric pressure. When 
the momentum of the Halon attempts to carry the Halon to a larger flow area, the 
pressure in the Halon jet starts to drop below atmospheric, and the surrounding 
air pushes the flow inward. 


■i 
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Multicone Nozzle 





Wide angle dispersion can be achieved only with discrete jets or sheets 
that penetrate the atmosphere but do not coalesce into a single flow stream. 
There must be air between the jets. The flow from a nozzle with hemispherical 
distribution will have fixed Halon flow area and increasing air-gap area as the 
Halon travels away from the nozzle. 



The problem with the showerhead nozzle was thought to be the straight 
orifices which discharged the Halon at throat pressure and allowed the Halon 
jets to spread without restraint and merge to exclude air gaps. To avoid this 
effect, a new nozzle was designed with converging-diverging passages. The 
diverging sections would give controlled, narrow-angle expansion to atmospheric 
pressure and produce jets that would penetrate the air without merging. 

A nested annular nozzle configuration was chosen for ease of fabrication. 
The resulting "multicone" nozzle, Figs. Ml and 42, has eight converging- 
diverging annular nozzles at angles of 15 to 85 degrees from center. The total 
throat area is 600 mm 2 and the exit area is 9000 mm 2 , giving an expansion ratio 







Figure k2. Multicone nozzle: (a) Parts 






of 15. The nozzle exit occupies 75 percent of the hemispherical exit surface, 
allowing an Initial 25 percent air gap area; at a distance of 1.0 m from the 

nozzle the Halon jets occupy only 1.4 percent of a hemispherical surface, and 
the distance between jet centers is 175 mm. 

nozzle is held together by two hoops. These hoops give an undesired 
interruption of the flow; a center support post could be used instead. 

Figure 43(a) shows the multicone nozzle operating with water and nitrogen 

flow as viewed from between the support hoops, and Figure 43(b) shows the flow 

as viewed edge-on to a support hoop. The flow shows little tendency to pull 

a cylindrical jet. However, the support hoops produce a 25 degree gap 
in the flow. 

Figure 44 shows movie frames of Halon flow from the multicone nozzle. The 
Halon spreads at a large angle, and the nozzle achieves the desired 
hemispherical distribution except for the gaps due to the support hoops. 

Figure 45 is a movie frame of a fire suppression test conducted by Dr. 
Harry T. Johnson at NASA White Sands Test Facility. Two multicone nozzles, one 
on the lower left and one on the lower right, spray Halon 1301 away from the 
camera and toward a diesel fuel c ire ln „ slmulated tank creu conlpartment _ ^ 

gaps m the flow due to the support hoops are very pronounced, otherwise, the 

Halon is spreading satisfactorily over a wide angle. The fire was extinguished 
in 150 ms. 


D. Personnel Impact Hazards 


Personnel struck by discharging Halon could be injured either by the total 
force or by the high velocity and low temperature of the flow. 


0 









7 b 











;ure HU . Hal on flow test of the multicone nozzle at 









Total force measurements were made using the apparatus shown in Fig. 46. A 
target plate mounted on a load cell was placed downstream of the bottle valve. 
The maximum force occurred at the start of discharge when the flow rate was 
about 15 kg/s. The isentropic Halon velocity at 5.2 MPa bottle pressure is 220 
m/s. The theoretical peak force is thus 15x220 = 3300 N = 740 lbf. The peak 
forces measured were higher, about 4000 N. With a wide-dispersion nozzle a 
person standing more than about 0.5 m from the nozzle would be struck by only a 
fraction of the flow, and the force would be much less. 

Based on past experience with water-and-nitrogen nozzles (Ref. 10), 
impingement of two-phase flow on bare skin is uncomfortable but not injurious. 
To simulate Halon flow conditions with water and nitrogen a nozzle with a 20 mm 
exit diameter was operated at an exit velocity of 140 m/s. The flow rate was 
0.6 kg/s and the void fraction of the jet was 0.99. Placing a hand in the jet 
gave a sensation of a bundle of needles being pressed into the skin, but there 
were no injurious effects. It can be concluded that impingement of Halon flow 
would not be injurious from the standpoint of drop penetration. The low 
temperature would be an additional hazard. 


E. Summary 

Halon can be dispersed by simple nozzles with orifices, slots, or 
deflectors if nothing more is desired than a few jets or sheets of Halon. 
Uniformly filling a large volume with Halon flow is more difficult, because 
numerous closely-spaced Halon jets or sheets are required and there must be air 
gaps between the streams. A multicone nozzle design achieved the desired air 
penetration of multiple Halon sheets, and gave a hemispherical distribution 
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except for gaps due to support struts. A center-supported multicone nozzle 
design would eliminate the gaps. 

The impact force of Halon leaving a bottle was 4000 N, but the force with a 
dispersion nozzle would be less. The impact of water-and-nitrogen flow on skin 
at a velocity similar to Halon discharge conditions was found to be uncomfor- 


table but not injurious 


o 
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APPENDIX A. DERIVATION OF THEORETICAL MODEL 


A. Assumptions 

The two-phase mixture of Halon liquid, Halon vapor, e.’-j nitrogen that flows 
out of the bottle is assumed to be a homogeneous mixture as it flows through the 
discharge piping and nozzles. The mixture behaves as a single-phase 
compressible fluid with a density that is a function of temperature and 
pressure, and all parts of the mixture have the same velocity. Inside the 
bottle the fluid has two layers, with the liquid confined to the lower layer. 

The flow in the bottle and pipes is assumed to be adiabatic, because the 
bottle and pipe wall areas are too small to transfer an amount of heat that 
would vaporize enough additional liquid to significantly affect the flow. 

A key simplifying assumption is that the discharge flow rate out of the 
system at each instant of time is equal to the steady-state flow rate that would 
exist if the bottle conditions were held constant at the values for that 
instant. The flow rate out of the bottle and into the piping system is then 
found by adding the rate of storage change in the piping system to the discharge 
flow rate, but the discharge flow rate itself is not a true transient 
calculation. This "quasi-steady" approximation is accurate during most of the 
discharge time, for piping volumes that are not greatly larger than the bottle 
volume. The approximation breaks down during initial pressurization of the 
piping system, when the storage rate dominates, and this process is treated by a 
different procedure. 
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B. Two-Phase Mixture Properties 

The Halon-nitrogen mixture is a two-phase, two-component mixture. The 
nitrogen will be designated "component A" and the Halon "component B», with 
subscripts a and b, respectively* 

By Henry’s Law, the mole fraction of nitrogen in the liquid phase is 

proportional to the partial pressure of tne nitrogen, p a . 

Thus, 

m a / w a p 

m aP 7W ae +( > - V )T V = HPa (1) 


where is the liquid mass, m a? is the mass of dissolved component A, W a£ and 
W bP are the molecular weights of A and B, respectively, and H is the mole 
fraction dissolved per unit pressure of A (the inverse of Henry’s Law constant). 
Solving for m a{; /m ? = a (the concentration of dissolved nitrogen) the result 


(W a£ /V W Hp a 
1 + ( ~ W a£ Tw b£- D'HPa 


The partial pressure p b 0 f component B is equal to the vapor pressure p bo of 
pure B multiplied by the mole fraction of B in the liquid (Raoult’s Law). 


P b = (1 - Hp a )p b0 


The total pressure is 


P = P a + P b 


A 




Combining Eqs. (3) and (4) the nitrogen partial pressure is 


P - P 


bo 


a ” 1 - Hp 


(5) 


bo 


The effective molecular weight W g 0 f the gas phase, such that the gas 
density is given by Wg p/RT where R is the universal gas constant, is 


W = 

e 


W p + W p 
ag a bg K b 


( 6 ) 


The mass fraction of Halon vapor in the gas phase is 


8 = 5 ^ 


(7) 


V 


The mass ratio of component A (gaseous nitrogen plus dissolved nitrogen) to 
component B (Halon liquid plus vapor) is the "component ratio" r c , defined by 


m , +• m 

r = _£i jag 

c m 0 + m. 
b£ bg 


( 8 ) 


The mass ratio of gas to liquid, in terms of r c and the concentrations a 
and [•. , is 


r - »(1+r ) 
c c 


m 1 - 6(1 +r ) 
c 


(9) 


0 1 


If P a } and are the densities of dissolved nitrogen and liquid Halon, 
respectively, then the density of the Halon-nitrogen solution is 



The specific enthalpies of the gas and liquid phases are the sums of the 


contributions from the components: 

h = (1-3) h + 3h. 
g ag bg 


arH 


h^ = ah g . + (1-c0h bJ 


( 12 ) 


The specific internal energies of the gas and liquid phases are found by 
subtracting the pressure-volume product from the enthalpy: 


u = h - 
g g 


(13) 


u, = hf 


(14) 


The two— phase mixture density, in terms of the phase densities, is 

1 * r_ 


mix 


m + 

3 g 


1 

p t 


(15) 


As a gi/en mass n t of the two-phase mixture proceeds through the flow 
system, the component ratio r c remains constant but the mass ratio r m changes 


and the gas and liquid flow rates change according to 

r 


m = 


m 


and 


_ ii«. 

g 1 ♦ r _ 

m 


m v = m t “ m g 


(16) 
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(17) 


C. Fluid Expansion in the Bottle 

At some intermediate pressure pi during discharge the mass of fluid in the 
bottle is mi and the volume of fluid is the bottle volume v bot - When the bottle 
pressure drops to a lower value P 2 the mass m-j expands to a larger volume equal 
to m 1 /P mix2 . where P mix2 is the density of the mixture (consisting of the liquid 
layer and the ullage gas) at pressure p 2 and temperature T 2 , with T 2 to be 
determined. The volume of fluid that must leave the bottle during the 
expansion from p-j to P2 is 

m 1 _ v (18) 
v out ” hot 

H mix2 

If the average density of the fluid in the liquid layer during the 
expansion from p-| to p 2 is o QUt , then the mass of fluid discharged is m 0 ut 
°out v out* To satisfy conservation of energy the change in the internal energy 
of the fluid in the bottle must equal the enthalpy of the fluid leaving. Thus, 


where h QU ^ is the specific enthalpy of the fluid leaving. 

The new temperature T 2 is the value that causes Eq. (19) to be satisfied. 
Once T 2 is found, all of the bottle conditions at the new pressure level p 2 are 
known, and the amount of fluid m out discharged to reach p 2 is known. 

Assumptions must be made about the location of the nitrogen bubbles in 
calculating the density c out and enthalpy h out of the fluid leaving the bottle. 
To decide the best assumption, predictions for Test 146 were made for two 
limiting cases: (1) no bubble rise, with all of the released nitrogen remaining 


trapped in the liquid layer, and (2) complete bubble rise, with all of the 
released nitrogen moving to the ullage. Fig. A-1 shows the comparison with data 
for the two assumptions. If complete bubble rise is assumed, the predicted 
discharge time is too short. The assumption of no bubble rise gives good 
agreement with the observed discharge time and is the assumption adopted in the 
model . 

In the model the dissolved nitrogen is held in solution until the nitrogen 
release pressure is reached. Fig. A-2, step 1-2, shows an example of the bottle 
conditions that satisfy Eqs. (18) and (19) for expansion without nitrogen 
release. If the nitrogen is then released from solution with no mass 
discharged, the conditions determined by Eq. (19) are shown in step 2-3. The 
next expansion step, with equilibrium dissolved nitrogen concentration, is shown 
in step 3-^. These illustrations show the detailed bookkeeping necessary in 
calculating two-phase, two-component expansion processes. 


D. Nozzle Flow 

In a nozzle the pressure drop is balanced by velocity increase. For a 
small pressure change dp the change in velocity is given by 


d(V 2 ) = 


2 dp 
Onix 


( 20 ) 


If Eq. (20) is integrated over a finite pressure decrease from pi to P 2 , 
and the mixture density is evaluated from Eq. 05) with gas density calculated 
from the perfect gas equation, a convenient equation for velocity increase, 
accurate for substantial pressure differences, is 


V 


2 

2 



1 ♦ r 


r p 
mm m 


gm 



Pl-P 2 


(21 ) 
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Figure A-l. 


Comparison, for no bubble rise and complete bubble rise, between 
experimental and theoretical pressure— time curves 
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Figure A-2. Typical bottle conditions at four successive stages of Halon 
discharge 


where r mm , p m , P gm , and P^ m are the arithmetic mean values of gas/liquid ratio, 
pressure, gas density, and liquid density, respectively, over the interval from 
Pi to P 2 . 

The temperature T 2 at pressure p 2 is found by guessing T 2 , using that guess 
in evaluating the mean values in Eq. (21), and then seeing if the enthalpy 
change of the fluid is equal to the kinetic energy change (V 2 ^ - V^)/2. 
Temperature T 2 is iterated until the kinetic energy change from Eq. (21) agrees 
with the enthalpy change. 

At the end of each pressure step the flow area for given flow rate m is 
found from 


O V 
*mix 2 2 


The smallest value of area found is the throat area A t . if the actual 
nozzle has throat area A, then the choked flow rate is mA/Af 


E. Pipe Flow 

In a pipe the pressure drop is balanced mainly by wall friction. For a 
pressure decrease from pi to P2 the velocity changes from Vi to V 2 , with V 2 
given by continuity: 


lmiX2 pipe 


The mean wall shear for the section of pipe between p-j and P2 is 


T v= 0,5 D mixm V m C f 


where ° m i xrn is the mean mixture density, V m * is the mean-square velocity, and Cf 
is the friction coefficient. 


U 


0 


The friction retarding force on the fluid in the length of pipe Ax between 
Pi and P2 is t w *Dpipe Ax. This force is balanced by the pressure drop and 
fluid momentum change. The pipe length needed for the required friction force 
is 


Ax = 


A . 
P*P e 


(p 1 - p 2 ) + m(V 1 -V 2 ) 

~T IT D~ 

w pipe 


(25) 


The new temperature T2 is found in the same way as for a nozzle, iterating 
on temperature until the enthalpy decrease equals the kinetic energy increase. 

As the flow proceeds down a long pipe the distance increment Ax required to 
provide enough friction to balance a given pressure difference pi - P2 
decreases. When a zero or negative Ax is found this means that sonic conditions 
have been reached. The pipe length at that point is the maximum length possible 
for the given flow rate. 


F. Abrupt Enlargement 

Between the valve and the nozzle or pipe inlet the flow area abruptly 
increases. The kinetic energy decrease downstream of the valve throat is 
converted to heat. To calculate the slight temperature rise downstream of the 
valve the downstream temperature is iterated until the kinetic energy decrease 
equals the enthalpy increase. 

G. Flow Rate Calculation for Series Components 

During discharge the Halon flows subsonically through several restrictions 
before reaching sonic conditions at the final restriction. For example, the 
Halon flows through the bottle valve and may then flow through a pipe inlet 
nozzle and a pipe before reaching the discharge nozzle throat. To calculate the 
flow rate through such a series string a flow rate must be guessed. The sonic 
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flow area required after the flow has encountered the various pressure losses 
can then be calculated. If the required area is larger than the actual throat 
area of the discharge nozzle then a smaller flow rate must be guessed. The flow 
rate is iterated until the required and actual throat areas agree. 

The final flow restriction moves from one component to the next during 
filling of the piping system. The flow-rate search is first done with the 
bottle valve area as the throat area, then with the pipe inlet orifice, then 
with successive positions in the pipe, and finally with the discharge nozzle. 
The position of the fluid front, which determines which flow restriction is 
controlling, is calculated from the fluid front (sonic) velocity and elapsed 
time. 


H. Elapsed Time Calculation 

For each pressure increment p-| to P 2 the average flow rate m^g at the 
discharge restriction is calculated as described above. The change in stored 
mass Am S £ 0 is found at the same time by integrating the fluid density over the 
pipe volume at the two pressures and finding the difference. The mass of fluid 
m out that left the bottle between pi and P2 is already known from the bottle 
calculations. The time increment At between pressures pi and P2 must be such 
that the sum of the discharged mass m<ji S At plus the storage increase Am s t 0 is 
equal to the mass leaving the bottle. Thus, the time increment between pi and 


P 2 is 


m . - '\m . 

out sto 


I. Pipe Pressurization 

Between the time the Halon reaches the end of the pipe and the time the 
pipe reaches peak pressure the storage change predominates and a different 
procedure for calculating elapsed time is used. The flow rate into the pipe at 


vi) 


the time the liquid reaches the end of the pipe is mini* ^he first flow rate 
out of the discharge nozzle, with just the pipe sonic exit pressure applied to 

the nozzle, is m ou t1* The mass stored in the pipe at this time is m s to1* 

The pressure at the time of peak pipe pressure is guessed and the 
corresponding discharge flow rate m out2 is calculated. Since the pipe pressure 
is not changing at that instant, having just reached its peak, the pipe inlet 

flow rate is equal to the outlet flow, m in2 = m 0 ut2* The mass stored at that 
time is m sto2 . 

The average inlet flow rate during pressurization is mj nm = (mj n l+ mi n2 )/2 
and the average outlet flow rate is m outm = (m 0 ut1 + m out2^ /2. 

The average net flow rate into the pipe is m^ nm _ m ou tm> an d that flow rate 


times the filling time Atfj .11 is equal to the storage increase m s to2 - m sto1* 
Thus, the time increment between the arrival of the liquid front at the end of 

the pipe and the instant of peak pipe pressure is 

m sto2 ~ m sto1 (27) 

'fill = . 

m . — m 

inm outm 

The mass Am bot discharged from the bottle during the time increment Atfin 
is the product of the average inlet flow rate minm and Atfin. From the bottle 
calculations the bottle pressure existing when that much mass has left the 
bottle is known and can be compared with the guessed bottle pressure at the time 
of peak pipe pressure. If the two pressures do not agree a new bottle pressure 
is guessed and calculations starting with ra out2 are repeated until the results 
converge. 


J. Multibranch Calculations 

For multibranch calculations ti ressure at the junction between the main 
pipe and the branch pipes is calculated from the previous flow rate. Each 
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branch is then calculated as though the junction conditions were the bottle 
conditions for that branch. The branch flow rates and branch storage masses are 
added and applied to the exit of the main branch and the calculations then 
completed in the same way as for a 3ingle-branch system. 

K. Summary 

The model consists of two sectionsi a bottle section that calculates bottle 
outlet conditions as a function of mass discharged and a flow section that 
calculates flow rate as a function of bottle outlet conditions. The flow rate 
calculations assume steady flow and determine flow rates from sonic flow 
conditions at the final flow restriction prevailing at the time. Adding the 
pipe storage changes gives the bottle discharge mass for each pressure 
increment, from which the elapsed time can be determined. For the pipe 
pressurization step the elapsed time is obtained from the average of the net 
flow rates at the beginning and end of pressurization. 



APPENDIX B. HALON AND NITROGEN PROPERTIES 



Thermodynamic properties of Halon 1301 are tabulated in Ref. 11 and 
Physical and transport properties in Ref. 12. Measurements of specific volumes 
and vapor pressures of solutions of nitrogen in Halon 1301 are reported in Ref. 
13. Based on Ref. 13, Du Pont developed a computer program called "SUPER" for 
calculation of nitrogen-pressurized Halon 1301 properties; the equations used in 
program SUPER are given in Ref. 14. 

Using the equations and tabulations of Refs. 11 to 14 and other sources, 
property subroutines were written for program HFLOW and for the modified version 
of SOLA-LOOP (Appendix C). The equations used are the following: 


1. 


Nitrogen Specific Volume 

The equation of state for gaseous nitrogen from Ref. 14 is 


P = , R » T - 99.164 

W(!v-0.015313) v(v+0.01$313VT 


( 1 ) 


where p is the nitrogen pressure in psia, T is the temperature in R, v is the 
specific volume in ft3 /lb , „ o ls th# unlversal gas conatant ^ ^ ^ 

ft3/mole-OR), and W is the molecular Height (28.016). The gas constant for 


nitrogen is R = R q /w = 0.38305. 


The equation of state can be rearranged for iterative solution as follows: 


AP = n-L 


v(v+ c i2)/T 


( 2 ) 


v = 


RT 


AP 


(3) 


where Cn = 99. 164 and C 12 = 0.015313. 


B-l 


0 1 


t> 


To start the solution, an initial guess at specific volume is calculated 


from the perfect gas law. 


v 0 = RT/p 


Then v Q is substituted in Eq, (2) to find the small pressure correction Ap f and 
a new volume v is found from Eq. (3). The answer converges rapidly with 
successive iterations. 

It is convenient to express the specific volume in terms of an effective 
molecular weight that gives the correct specific volume when substituted in the 
perfect gas equation. The effective molecular weight is 


» — wy 

6 vp 

2. Dissolved Nitrogen Specific Volume 

The partial specific volume of nitrogen dissolved in Halon 1301 is given by 


Ref. 14 as 


v = C 21 + C 22 2-n (C 23 - t) 


where v is the increase in volume of a nitrogen-Halon solution in ft 3 per lb of 
dissolved nitrogen and t is the temperature in °F. The constants are C 21 = 
0.24094456, C 22 = -0.04021084 and C 23 = 152.96. 


3. Nitrogen Enthalpy 

The enthalpy of nitrogen over the range 0 to 600 psi and -80° F to 80° F 
can be represented by 

h g = C 31 P ♦ C 32 (300-T) + C 33 p(300-T) + C 34 p(300-T) 2 (7) 

where h g i s the enthalpy in J/g-mole, p is the pressure in at',, and T is the 
temperature in K. The constants are C 31 = -5.8462, C 32 = -28.96, C 33 = -0.0428, 
and C 34 = -0.3788 x 10“ 3 . 





4. 


Halon 1301 Saturation Pressure 




\\ 


The saturation pressure of Halon 1301 from Ref. 14 is 

loglC) P = Cn + C 12 /T 2 + C 13 /T + C 14 T + C 15 T 2 


( 8 ) 


where p is the pressure in psi and T is the temperature in R. The constants are 

C 11 2 7.322506, C 12 = -44295.0, C 13 = -1795.678, C 14 = -0.00438339, and C 15 = 
0.312584 x 10-5. 


5* Halon 1301 Vapor Density 

From the Ref* 11 tabulation of saturated vapor density an effective 
molecular weight can be calculated that varies from 153.77 at -76°F to 202.02 at 
76°F (the actual molecular weight of Halon 1301 is 148.93). When the effective 
molecular weights are substituted in the perfect gas law together with the 
saturation pressure from Eq. (8), the density agrees with the Ref. 11 tabula- 
tion. A curve fit to the effective molecular weight that agrees within 0.1 
percent is given by 

W g = C 2 i + C 22 t + C 23 t 2 + C 24 t 3 (9) 

where t is the temperature in °C. The constants are C 2 i = 178.10867, C 22 = 
0.74214272, C 23 = 0.008111179, and C 24 = 0.42822036 x 10" 4 . 


6 . Halon 1301 Liquid Density 

From Ref. 14 the saturated liquid density is 

P i = c 31 * C 32 (T c -T) + C 33 (T c -T) ^ 2 + C 34 (T c -T) ^3 + c 35(Tc _T)2 Qo) 

where is the density in lb/ft 3 , T is the temperature in R and T c is the 

critical temperature of 612.27 R. The constants are C 31 = 46.50, C 32 = 
0.039808, C 33 = 0.82202, C 34 = 9.30627, and C 35 = 0.922 x 10 -6 . 
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Halon 1301 Vapor Enthalpy 

A curve fit to the Ref. 11 tabulation of saturated 
agrees within 0.04 percent is given by 


vapor enthalpy that 


h g = C 41 + c 42 fc ♦ c 43 fc2 + C 44 t3 ( 11 ) 
where h g i 3 the enthalpy in Btu/lb-°F and t is the temperature in °C. The 
constants are C 41 = 53.635704, C 42 = 0.10834558, 0143 = -0.69031717 x 10 - 3 , and 
c 44 = -4.1710027 x 10 -6 . 


® • Halon 1301 Liquid Enthalpy 

A curve fit to the Ref. 11 tabulation of saturated liquid enthalpy that 
agrees within 0.1 percent is given by 

h l = C 51 + C 52 t + C 53 t2 + c 54 t3 (12) 

where h £ is the enthalpy in Btu/lb-°F and t is the temperature in °C. The 

constants are C 51 = 12.690278, C 52 = 0.33714354, C 53 = 0.5673509 x 10"3 and c 5l4 
= 1.671468 x 10 - 6 . 


9. Henry's Law Constant 

From Ref. 14 the Henry's Law Constant is given by 


H = C 61 ♦ Cg 2 t + Cg 3 t 2 + c 64t^ 


(13; 


where H is the partial pressure of nitrogen (in psia) per mole fraction of 
dissolved nitrogen, and t is the temperature in °F. The constants are c 61 
4763.8683. C .2 = 3.9*164387, C b} = -0.016242350, and c 6fl = -0.4445893 x 10-3. 


10. Halon 1301 Liquid Viscosity 


The Ref. 12 tabulation of Halon 1301 liquid viscosity can be represented 


by 


= C.112 + 3. 11x10“5( 150-t) 


1.673 


( 14 ) 


where U£ is in centipoise and t is the temperature in °F. 


11. Halon 1301 Vapor Viscosity 

The Ref. 12 tabulation of Halon 1301 vapor viscosity can be represented 
by 


U g = 0.0137 + 3x10 _5 t 


( 15 ) 


where t is the temperature in °F. 


12. Halon 1301 Surface Tension 

An equation that follows the fourth-power rule and fits the Ref. 12 data 
at 0°F and 70 °F is 

o= 6.6x10“8(p£-Pg)4 ^ 

where a is the surface tension in dyne/cm and o, and P g are the liquid and vapor 
densities in lb/ft^. 


13. Conversion to SI Units 

The conversion factors to SI units for these properties are: 
psi x 6894.76 = Pa 
lb/ft3 x 16.01846 = kg/m^ 

Btu/lb x 2326.0 = J/kg 
ce.itipoise x 10~3 = N-s/m^ 
dyne/cm x 10“3 = N/ m 
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APPENDIX C. FLOW CALCULATIONS WITH THE SOLA-LOOP PROGRAM 


A group at the Los Alamos National Laboratory (LAND developed a computer 
program called SOLA-LOOP for calculation of transient two-phase flow in piping 
networks (Ref. 15). The program was written for steam and water. We prepared 
a modified version of SOLA-LOOP with Halon 1301 and nitrogen properties 
replacing the steam and water properties. The modified program was run on the 
LANL computer from a terminal at JPL. 

The SOLA-LOOP program divides the piping system into many small segments 
and solves the continuity, momentum, and energy equations for each segment at 
each time step. There are no restrictions on the geometry of the piping system. 
With sufficiently small segments and small time steps the results should be 
highly accurate. 

Predictions of the modified SOLA-LOOP program were compared with test 
results for the single-branch system of Fig, 18. Figure C-1 shows the results 
for Test 174, with a 77 mm^ nozzle. The agreement is very good. The SOLA-LOOP 
program not only follows the average pressure aecay curves accurately but shows 
the pressure oscillations following the filling of the pipe, although with some 
exaggeration. 

Figure C-2 compares the theoretical and experimental pressures for Test 
177, with no orifice on the end of the pipe. The SOLA-LOOP program correctly 
shows the gradual pressure rise at the pipe exit, because the model includes the 
effect of air in the pipe ahead of the flow. 
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TEST 174 
BOTTLE PRESSURE 
77 mm^ NOZZLE 

MEASURED 

THEORETICAL 


TIME, s 


Comparison, for discharge through a pipe and nozzle, between 
(a) Bottle pressure SOLA-LOOP predictions and data 






PRESSURE, MPa 








Computing time and cost for the SOLA-LOOP program on the LANL CDC-7600 
computer ranged from 30 s and $4, for computing the flow from a bottle through 
an orifice, to 11 minutes and $90 for a three-branch distribution system. 

An improved version of the SOLA-LOOP program, called "SOLA-NET", is now 
available (Ref. 16). Adding Halon properties to that program would probably 
give the most accurate Halon flow program possible at the present. 


APPENDIX D. TEST SUMMARY 


In addition to the tests run for comparison with the theoretical models, 
about 150 tests were made with flow through the bottle valve only. These tests 
were referred to as the "short system tests". Two different bottle sizes were 
used, with valves from four different manufacturers. Four different Halon fill 
densities were used. Table D— 1 lists the short-system test numbers for each 
combination of conditions. 

The bottle pressure used for each short-system test was the value that 
would have been reached on heating or cooling a standard-condition bottle (70°F 
and 750 psig) to the test-cell temperature. That pressure was obtained from the 
graphs shown in Figs. D-1 (a) and (b). 

The results of each test consisted of plots of pressure, temperature, and 
target force versus time. These plots are voluminous and are available 
separately. 

A summary of test conditions and results is presented in Table D-2 for the 
short-system tests and also for the tests used in developing the models. An 
estimate of liquid discharge time is given for those tests where a slope break 


could be identified 


PRESSURE, psia 
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D-l. Curves used for setting bottle pressures in the short-system test 
(a) With a 170-in^ bottle 
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Table D-1. Short-System Test Numbers 


HRM 

VALVE 


MAROTTA 

VALVE 
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Test Summary 


a 



Test 

No. 

Valve 

Bottle 
Volume , 
in. 3 

Halon 

Fill, 

lbm 

Fill 
Density, 
lbm/ ft 3 

initial 
Pressure , 
psia 

Initial 
Liquid 
Temp ,°F 

Nitrogen 
Release 
Pressure , 
psia 

Liquid 
Expulsion 
Time , 
ms 

46 

CR-1 

169.0 

4.98 

50.9 

785 

85 

605 

90 

47 

CR-1 

169.0 

5.02 

51.3 

7^3 

65 

460 

100 

49 

CR-1 

169.0 

5.25 

53.7 

750 

62 

432 

100 

50 

CR-2 

169.0 

5.37 

54.9 

787 

75 

570 

— 

51 

CR-2 

169.0 

5.05 

51.7 

771 

72 

700 

— 

52 

CR-2 

169.0 

6.15 

62.9 

725 

6l 

455 

— 

53 

CR-2 

169.0 

5.9^ 

60.7 

751 

64 

460 

150 

54 

CR-2 

169.0 

5.91 

60.4 

752 

66 

480 

— 

55 

CR-2 

169.0 

6.08 

62.2 

748 

68 

538 

— 

56 

HR-1 

231.5 

6.98 

52.1 

769 

68 

525 

150 

57 

HR-1 

231.5 

7.04 

52.5 

747 

72 

490 

l4o 

58 

HR-1 

231.5 

7.07 

52.8 

782 

74 

541 

160 

59 

HR-1 

231.5 

7.13 

53.2 

769 

74 

549 

160 

60 

HR-1 

231.5 

5.06 

37.8 

775 

78 

561 

100 

6i 

HR-1 

231.5 

5.02 

37-5 

778 

75 

559 

110 

62 

HR-1 

231.5 

5.05 

37-7 

745 

67 

490 

— 

63 

rIR-1 

231.5 

5.08 

37-9 

763 

67 

525 

— 

64 

HR-2 

231.5 

6.99 

52.2 

737 

64 

464 

150 

65 

HR-2 

231.5 

T .00 

52.3 

739 

65 

465 

160 

66 

HR-2 

231.5 

7.00 

52.3 

733 

66 

480 

150 

67 

HR-2 

231.5 

6.98 

52.1 

750 

65 

490 

150 

68 

CR-1 

229 . 3 

5.08 

37-9 

717 

60 

430 

80 

69 

CR-1 

229.3 

5.05 

37-7 

744 

64 

460 

80 

70 

CR-1 

229.3 

4.93 

36.8 

720 

63 

4oo 

80 

71 

CR-1 

229.3 

5.07 

37.8 

732 

61 

4ii 

80 

72 

CR-1 

229.3 

6.98 

52.6 

739 

63 

620 

150 

73 

CR-1 

229.3 

7.10 

53.5 

736 

62 

451 

150 

74 

CR-1 

229.3 

7.10 

53.5 

733 

62 

450 

150 

75 

CR-2 

229.3 

7.02 

52.9 

650 

53 

265 

— 

76 

CR-2 

229 . 3 

6.88 

51.8 

737 

63 

475 

— 

77 

CR-2 

229.3 

6.88 

51.8 

737 

62 

450 

— 

78 

CR-2 

229 . 3 

7.05 

53.1 

781 

67 

535 

150 

79 

CR-2 

229.3 

6.97 

52.5 

735 

67 

483 

iko 

81 

HT-1 

229.6 

5.04 

37.9 

732 

59 

433 

90 

82 

HT-l 

229.6 

5.05 

38.0 

761 

66 

491 

80 

83 

H.T-1 

229.6 

5. Ll 

38 . 5 

745 

64 

•V70 

80 

84 

HT-1 

229.6 

5.LS 

39.0 

751 

64 

4 1 4 

8u 

85 

CR-2 

229 . 3 

6.92 

52. L 

968 

113 

— 

— 

85,, 

CR-2 

229 . 3 

7.07 

53.2 

1043 

1 1 f 

970 

— 

87 

CR-2 

229 . 3 

7.17 

54 .0 

J.017 

123 

911 

— 

88 

CR-2 

229 . 3 

6.99 

52 . 4 , 

1022 

125 

900 

— 

89 

CR-2 

229.3 

7 . 00 

62.7 

1005 

1.24 

9 bh 

— 

90 

CR-1 

229 . 3 

7.0L 

52.8 

985 

119 

9^1 

— 

91 

CR- L 

229 . 3 

7.13 

53.7 

1003 

1 24 

909 
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Table D-2 . Test Summary (Contd.) 


Test 

No. 

Valve 

Bottle 
Volume , 
in. ^ 

Halon 

Fill, 

Ibm 

Fill 

Density , 
lbm/ft.3 

Initial 
Pressure , 
psia 

Initial 
Liquid 
Temp ,°F 

Nitrogen 
Release 
Pressure , 
psia 

Liquid 
Expulsion 
Time , 
ms 

92 

CR-l 

229.3 

7.13 

53.7 

982 

123 

— 

— 

93 

CR-1 

229.3 

7.15 

53.8 

1028 

120 

980 

— 

9k 

HR-2 

231.5 

7.01 

52.3 

989 

121 

— 

— 

95 

HR-2 

231.5 

6.97 

52.0 

1034 

123 

945 

— 

96 

HR-2 

231.5 

7.03 

52.5 

969 

120 

— 

— 

97 

HR-2 

231.5 

7.10 

53.0 

1050 

120 

974 

— 

98 

HR-1 

231.5 

6.81 

50.8 

1007 

124 

958 

— 

99 

HR-1 

231.5 

7.04 

52.5 

999 

119 

— 

— 

100 

HR-1 

231.5 

7.05 

52.6 

975 

120 

950 

— 

101 

HR-1 

231.5 

7-08 

52 . 8 

997 

122 

970 

— 

102 

HR-1 

231.5 

6.92 

51.7 

802 

65 

581 

170 

103 

HR-1 

231.7 

6.97 

52,0 

781 

67 

560 

— 

104 

HR-1 

231.7 

7.18 

53.5 

797 

76 

600 

— 

105 

HR-1 

231.7 

7.14 

53.5 

848 

70 

682 

190 

106 

HR-1 

231.7 

7.07 

52.7 

841 

80 

680 

180 

107 

HR-1 

170.7 

5.01 

50 . T 

805 

81 

634 

0 

0 

H 

108 

HR-1 

170.7 

5 . 01 

50.7 

839 

88 

687 

100 

109 

HR-1 

170.7 

4.99 

50.5 

796 

79 

624 

110 

110 

HR-1 

170.7 

4.99 

50.5 

801 

78 

630 

100 

111 

HR-1 

170.7 

5.10 

51.6 

775 

70 

550 

120 

118 

HT- 1 

168.2 

5.95 

6l . 1 

820 

83 

667 

150 

119 

HT-1 

168.2 

5.96 

61.2 

821 

84 

675 

l4o 

121 

HT-1 

168.2 

5.12 

52.6 

783 

74 

590 

90 

122 

HT-1 

168.2 

4.99 

51.3 

826 

85 

680 

80 

123 

HT-1 

168.2 

4.97 

51.1 

830 

89 

680 

80 

124 

mr-i 

168.7 

6.12 

62.7 

853 

91 

730 

200 

125 

MR-1 

168.7 

5.88 

60 . 2 

771 

73 

560 

190 

126 

MR-1 

168.7 

5 . 9 c 

60 . 4 

790 

77 

605 

190 

1.27 

MR-1 

168.7 

5.91 

60 . 5 

787 

81. 

570 

180 

128 

MR-1 

168.7 

4.98 

51.0 

806 

81 

631 

iho 

129 

MR-1 

168.7 

4.99 

51 .1 

802 

80 

632 

l4o 

130 

MIR- 1 

168 . 7 

4 . 9 3 

50.5 

806 

77 

6 3.1. 

140 

.1.3.1 

MR- L 

168.7 

9.0 2 

5 l.Y 

751 

72 

‘♦79 

1-30 

L32 

i LI-.— 1 

170.7 

9.92 

60.0 

780 

76 

5 \ 'j 

MfO 

183 

HH- 1. 

170.7 

9 . 96 

6 0 . 3 

794 

67 

586 

L7-.' 

i 84 

HR— ] 

170.7 

# )■) 

9 9 • i 

83! 

8 1 j 

654 

l.yd. 

L 

HR-1 

l7o.7 

v : . 9 i 

i i 1 j , 1 ; 

793 

74 

6 1.0 

lv ' 

!_ ^ t , 

: i ; - 1. 

JmB..! 

9 .'R-o 

r o.O 

790 

77 

6 L 5 

.1.4' ) 

LTf 

MR-F 

8 0 .;. 

0 . 9 r 

: > 1 

1 

823 
y. ^ , 

80 

7 

1 >80 

( . ' "i 

IF. 1 » 

].3b 

H : - 1 

j ; ■ ■ i 

) * f > 
FFR.o 

T.r-J 

y - • -I 
9 r 7 . 4 

4 1 

- j . 1 . 

1 ✓ 
78 

* V ‘ 

F'.jO 

lOo 

HI- L 

220.6 

7 . 39 

> > • L - 

O' 3 9 

78 


.1 40 

L\] 

H ; 

220 . 6 

7.'^ 

53.1 

8 1 b 

82 

o3Y 

1> j 

I.-O 
1 V- 

Hl-i. 

i » i'*i , y 

4 . 36 

44 .8 

i 

83 

019 

— 

H 1-1 

1 . 2 

9.13 

{ ' ') " / 
y 1 - • i 

75'' 

71 

9 3 n 

8:0 
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T e s t Summary ( C on t d . ) 


Table D- 









Nitrogen 

.Liquid 



Bottle 

Halon 

Fill 

Initial 

Initial 

Release 

Expulsion 

Test 


Volume , 

Fill, 

Density , 

.Pressure , 

Din aid 

Pressure 

. Time « 

No. 

Valve 

in. 3 

.lbm 

lb.'i./ft , 3 

psia 

!e:;ip ,°i 

psia 

ms 

150 

HT-2 

229.6 

7.06 

53.1 

813 

77 

635 

160 

151 

HT-2 

229.6 

6.94 

52.2 

811 

77 

625 

160 

152 

HT-2 

229.6 

7 . 02 

52.8 

806 

72 

629 


153 

HT-2 

229.6 

7.01 

52.8 

970 

155 



154 

PIT-2 

229.6 

6.99 

52.6 

940 

108 

800 

]6o 

155 

HT-2 

229.6 

6.82 

5 i . 3 

800 

101 

670 

.1 60 

156 

HT-2 

229.6 

7.16 

> d . 9 

962 

124 




157 

HT-2 

229 . 6 

7 . 00 

52.7 

990 

160 



158 

HT-2 

229.6 

6.97 

52.5 

975 

120 




159 

HT-1 

168.2 

5.94 

61.0 

880 

128 

790 


160 

PIT-1 

168.2 

5.03 

5.17 

880 

105 

775 

90 

166 

MR-2 

229.4 

7.15 

53.9 

764 

71 

561 


167 

HT-1 

168.2 

5.17 

53.1 

777 

71 

571 

100 

168 

HT-1 

229.6 

7 . 01 

52.8 

790 

74 

605 

150 

169 

HT-1 

229.6 

7 . 02 

52.8 

1012 

125 

950 


170 

HT-1 

229.6 

7.01 

52.8 

1021 

159 

961 


171 

PIT-1 

229.6 

6.98 

52.5 

970 

125 

905 


172 

HT-1 

229.6 

7.02 

52.8 

985 

129 

930 



173 

HT-1 

229.3 

7.03 

53.0 

764 

68 

550 

800 

17U 

HT-1 

229.3 

6.95 

52.4 

76l 

80 

588 

2000 

175 

HT-1 

229.3 

7.01 

52.8 

756 

79 

590 

5d0 

177 

HT-1 

229.3 

7.02 

52.9 

757 

18 

586 

360 

i J 

MR-1 

229.4 

5.02 

37.8 

723 

54 

577 

100 

276 

MR-1 

229.4 

4.96 

37.4 

735 

61 

470 

90 

277 

MR-1 

229.4 

5.10 

38.4 

736 

61 

479 

90 

278 

MR-1 

229.4 

4.96 

37.4 

767 

6.1 

520 

100 

279 

MR-1 

229.4 

7.03 

53.0 

746 

62 

593 

200 

280 

MR-1 

229.4 

6.98 

52.6 

736 

65 

507 

200 

281 

MR-1 

229.4 

7.01 

52.8 

744 

65 

518 

200 

282 

MR-1 

229 . 4 

7.05 

53.1 

701 

61 

446 


283 

MR-1 

229.4 

7.04 

53.0 

1149 

145 

800 

200 

284 

MR-1 

229.4 

7.13 

r 

70 l 

1210 

151 

740 

170 

285 

MR-2 

229 . 4 

0.95 

52.4 

746 

65 

480 



286 

MR -2 

229 . 4 

6.95 

52.4 

754 

67 

510 



287 

MR-2 

229.4 

6.99 

2.4 

751 

69 

512 



288 

MR-2 

229.4 

7.04 

53.0 

745 

69 

510 

^ 

289 

MR-2 

229 . 4 

7.23 

4.5 

795 

81 

600 

210 

290 

MR -2 

229.4 

4.99 

37 • 6 

1490 

146 

1380 



298 

HT-1 

229.3 

7.00 

52.5 

783 

75 





294 

HT-1 

229.3 

7.05 

53..] 

761 

72 

— 

— 
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